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Adiabatic limit is the presumption of the adiabatic geometric quantum computation and of the 
adiabatic quantum algorithm. But in reality, the variation speed of the Hamiltonian is finite. Here 
we develop a general formulation of adiabatic quantum computing, which accurately describes the 
evolution of the quantum state in a perturbative way, in which the adiabatic limit is the zeroth-order 
approximation. As an application of this formulation, non-adiabatic correction or error is estimated 
for several physical implementations of the adiabatic geometric gates. A quantum computing process 
consisting of many adiabatic gate operations is considered, for which the total non-adiabatic error is 
found to be about the sum of those of all the gates. This is a useful constraint on the computational 
power. The formalism is also briefly applied to the adiabatic quantum algorithm. 



Recently a considerable amount of attention has been 
paid to the idea of using geometric phases accumulated 
by an adiabatically time-dependent Hamiltonian to re- 
alize quantum gates IM- Construction of universal 
gates by geometric quantum teleportation was studied, 
with the analysis of errors from imperfect control Q . On 
the other hand, adiabatic evolution is also the basis of the 
so-called adiabatic quantum algorithms , for which the 
speed and the overall time have been analyzed . 

For these quantum computing schemes to work, it was 
supposed that the adiabatic limit is retained. However, 
in practice, and particularly in the case of quantum com- 
putation, where the advantage lies in speedup and the 
operation time should be shorter than the decoherence 
time, the evolution is required to be completed in a fi- 
nite period of time. Therefore, it is important to know 
the full picture of the evolution of the quantum state 
and the non-adiabatic correction, which gives rise to er- 
ror if the adiabatic limit is necessary for the designed 
quantum computing scheme. Here we develop a general 
formulation of adiabatic quantum computing, applicable 
to the previously proposed quantum computing schemes 
and to any slowly varying Hamiltonian. As an adiabatic 
perturbation theory, it accurately describes the quantum 
evolution in a perturbative way, in which the adiabatic 
limit is the zeroth-order approximation. As an applica- 
tion, an examination is made on the non-adiabatic errors 
in several previously proposed implementations of the 
adiabatic geometric gates. We also investigate the non- 
adiabatic error in an entire quantum computing process 
consisting of many adiabatic gates, which has not been 
considered previously. Finally we briefly discuss the adia- 
batic quantum algorithm, noting that such an algorithm 
can still be implemented even if the non-adiabatic cor- 
rection is not vanishingly small. 

If the evolution of a time-dependent Hamiltonian is suf- 



ficiently slow, the adiabatic theorem tells that in the adi- 
abatic limit and under such conditions as continuity, non- 
crossing and differentiability, an instantaneous eigenstate 
at an initial time evolves to a state close to the corre- 
sponding instantaneous eigenstate at a later time 

In general, using the instantaneous eigenstate |0„(t)), 
one can always expand the state of the system \ip{t)) as 

IVW) -5]«"WI'^"W)e^PNnW], (1) 
n 

where ??n(i) = ^ji Jq En{T)dT is the dynamic phase. 
Then the Schrodinger equation ihdt\ip{t)) = H{t)\ijj) 
leads to 

dtan{t) = -^a„(t)(0„(i)|9i0„(t))exp[i?7„(t)-i?7„(t)], 
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(2) 

which, together with the initial condition a„(0) = 
((/)„(O)|'0(O)), determines \ip{t)). 

First suppose |V'(0)) is a non-degenerate eigenstate 
|0„(O)). Then in the adiabatic limit, one obtains 
IV'(0> ~ \(j)n{t)) exp[i^ri{t) + irinit)], where 7„(i) = 
/o('/'«|5r0«)dT = j{4>n\d^(t)n)dx'^ IS the geometric or 
Berry phase j^]. 

However, when non-adiabatic correction is considered, 
the exact state should be the solution of Eq. ||2Jl. For a 
slowly varying H(t), using a perturbative approach, one 
can obtain 



[/(t)|0„(O)) = exp[i7„(t) +i77„(t)][|0„(t)) + 



(3) 



S,„(t)-B„(t) 



In general, as in quantum computing, |^/'(0)) = 
a„(0)|(/)„(0)) is a superposition of different eigen- 
states. Then linearity of quantum evolution implies 



=^a„(O)C/(t)|0„(O)), 



(4) 
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where each U{t)\(j)n{0)) is as given in Eq. There- 
fore, an{t)exp[if]n{t)] = J2ni{^n{t)\U{t)\(f>^{0))am{0)- 

From ©, one obtains ((/>„(t)|t/(<)|0„(O)) « 
(0„(t)|t/(o)(t)|<^„(O)) = exp[i7„(i) + i7]nit)], while for 
n ^ m, (</.„(t)|C/(t)|0™(O)) ^ (0„(<)|C/(i)W|(/)™(O)) = 
?iexp[i7m(t) + ir]jn{t)]{(j)n\dt<j)m}/[En{t) - E„i{t)]. Here 
refer fcth order term. Since En{t) ^ Emit), 
Vn (t) 7^ Vm (t) ■ If one implements an all-geometric gate, 
in which the instantaneous basis states are \(j)n(t)) and 
\<t>mit)), the difference between rjnit) and Vmit) needs to 
cancelled out by using a certain method 

In the presence of degeneracy of eigenstates, denote the 
eigenstates as 100,^(0)), where n labels the energy levels, 
while an labels the different eigenstates in the subspace 
n. As the generalization of Eq. Q, we obtain 



C^WI«„(0))=exp[*,7„(i)][|xS„W>- 

tV- lxg'„(*))(xg'„(*)|3tX^„(t)) i , 

" l^m^n 2^I3,„ E^(t)-E„{t} \ "* 



(5) 



with ixs„w> = j:p„vzpjt)\rpM' ^"W = 

Pexp/„*A"(r)dr, where A^^^^ = (x^„(T)|5tx2„(^)> is 
the connection in the subspace n. V"' may be called non- 
abelian geometric phase or Wilczek-Zee (WZ) phase 0. 
In the zeroth order, U{t) is block-diagonal, each block 
being a WZ phase in the subspace of a set of degenerate 
eigenstates. In the adiabatic limit, as a unitary transfor- 
mation, a non-abelian geometric phase, i.e. the first term 
in ISJ, may be used to realize a quantum gate 

With the existence of degeneracy of eigenstates, a gen- 
eral superposition state can be written as 

IV'W) ^T.T.<MxUt))^M^VrM- (6) 

n an 

By choosing an appropriate basis for each degenerate 
subspace, the initial state can always be expanded in 
such a way that its projection in each degenerate sub- 
space is a single eigenstate I^Jj (0)), i.e. |V'(0)) — 
E„«^„(0)|x^„(0)), with |x^„(0)) = ir^jO))- Therefore 



(7) 



Thus it)exp[iijnit)] 

En^EpJx'^MU(t)\x'^JO))a'^JO)■ From ©, 

(xS„W|c^Wlx^„(o)) « {x'ijt)\u^"\t)\xxm) - 

e'xjp[ir]n{t)]Sa , while for n ^ m, 

(xS„Wlc^Wlx^„(o)) « {xUm^'Ht)\xZi^)) = 

hexp[iij„,mxZ,MdtxZit))/[Enit) - Emit)]. 

Through this formulation, it becomes clear that 
the adiabatic quantum computing is based on 
(Mm^'KmniO)) or (xS„W|t^(°Hi)lx^„(0)), with 
higher-order terms neglected. Besides, while the previ- 
ous proposals of adiabatic geometric gates are based on 
closed paths, there is nothing in principle against using 
open paths, as far as the corresponding geometric phases 
can be detected Another noteworthy point, which 



was not pointed out before, is that when the qubits 
under a gate operation is entangled with other qubits, 
the linearity of quantum evolution guarantees that the 
gate operation is still given by Eq. ((SJ or where the 
eigenstates are those of this concerned gate; one may 
include in the coefficients a„(0) or aJJ^(O) the states of 
the other qubits projected in the same branch as the 
eigenstates of the gated qubits. This is crucial for the 
possibility that different adiabatic geometric gates can 
be networked. 

There is a significant difference in the uses of (abelian) 
Berry phase and WZ phase to realize a quantum gate, 
under adiabatic limit. For a Berry phase gate, it is nec- 
essary to have d non-degenerate states, where d is the 
Hilbert space dimension of the gate. For a WZ phase 
gate, one intentionally restrict the gate in a single degen- 
erate eigenspace. A quantum gate based on WZ phase is 
more advantageous than that based on Berry phase, on 
the aspect that for the former, in the adiabatic limit, the 
state is always an instantaneous eigenstate of the Hamil- 
tonian, hence there is no dynamical phase difference be- 
tween the basis states, and it is more stable against en- 
vironmental perturbation. 

The non-adiabatic correction or error at time t is 

eit) ^ Pit) - u(°Hmm) = EZiU^'KtMm « 

;7(i)(t)|'0(O)). The adiabatic limit means |e(r)| < 1. 
For a Berry phase gate, with |-0(O)) = J2n^niO)\4'niO)), 

eit) « ^a„(0) J2 {'l^n.it)\U^'Ht)\KiO))\<l^niO))- (8) 
For a WZ phase gate geometric gate at En, with \ip{0)) = 

Ea„ aa„(0)|C„(0)), 



<t) 



a«„(0) 



EE 



(xLW|t^^'^WIxa„(o))IC„(o)). 

(9) 

Note that the first order correction at time t is de- 
termined only by eigenvalues, eigenstates and their time 
derivatives at t, hence is history- independent. This sim- 
plifies the analysis. The time derivatives do depend on 
the details of time-dependence. However, since only the 
path is specified , without the fine control of the dy- 
namics, numerically it suffices to obtain the order of mag- 
nitude. The first-order correction is ^ h/AT, where T 
is the time duration of the gate operation, A is the min- 
imum energy gap with other eigenstates. It is the pre- 
sumption of "slow variation" or perturbative approach 
that h/AT < 1. The A:-th order correction is ~ iU/AT)^. 

As applications, we now apply the above results to sev- 
eral physical implementations previously proposed. The 
first proposal, based on Berry phase, uses NMR Q. The 
Hamiltonian, in the rotating frame, is Hit) — R(t) • I, 
where R = iRx,Ry,Rz) — cos sin 0, cjq — ^), 
I = ^((Ta;, (Ty, (Tz), ujQ is proportional to the static mag- 
netic field in z direction, uji is proportional to the RF 
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magnetic field in xy plane, lo is its angular frequency, 
(j) is its initial phase. The instantaneous eigenstates is 

\l{t)) = l/VR[iR. - iRy)/VR^\ T) + VR^\ I 

)], with eigenvalue R/2, and \0{t)) ^ 1/VR[{R^ - 
iRy)/^/RTRl\ T)-/RTi?r| i)] with eigenvalue -R/2, 
where R = |R|. From this, one obtains, for n = 0,1, 
UnniO^t) = e^p[ijn{t) + irinit)]. r?o(i) = Rt/2, r]i{t) = 
—Rt/2. The Berry phase 7n(i) is, in the case of a cycle 
path C, n{C)/2 for n and -f](C) for n = 1, where 
r2(C) is the solid angel that C subtends at i? = 0. It is 
straightforward to write down Uoi{0,t) and C/io(0,t). For 
a gate operation of a period T, the order of magnitude 
of these two matrix elements, and thus the non-adiabatic 
correction, is h/RT = l/y/ojl + {uiq — uYT. The two- 
bit gate, of qubits a and b, is effected by addition of 
the interaction hJIazhz- For the conditional phases of 
qubit a, Wao shifted to to ujao + Jhz = Woo ± J/2, de- 
pending on the basis state \Ibz) of b. One can obtain the 
non-adiabatic corrections in the two subspaces, with the 
substitution of Waoi<^/2 for in R above. For a gate as 
in P, the gaps are of the order of several hundred Hertz, 
while T is of the order of second, hence the non-adiabatic 
corrections are of the order of 10~^. 

This method was also applied to a Josephson junction 
circuit The effective Hamiltonian is still as that 

for NMR, now with R = cos a, —iJj sin a, i?c(l — 
2noff)), where Ej and a are decided by the Josephson 
couplings of two junctions, is charging energy, 2enoff 
is the offset charge. In the charging regime, as used, 
Ec > Ej. Thus the non-adiabatic correction is of the 
order of h/^/Ej + E^{1 - 2noff)^T. Hence if 1 - 2noff 
is not too small, the adiabatic condition is Ti/EcT <C 1, 
more relaxed than previously thought 0| . 

An implementation of WZ phas e g ate was pro- 
posed for trapped ions 0, Q, Il5l |^ The one- 
bit gates are based on the Hamiltonian H = 
?i|e)(a;o(0| -I- wi(l| + i^a{a\ + h.c), One can find that the 
eigenstates are: = (ci;|e)-|-cjQ|0)-|-a;i]'|l) -|-a;*|a))/V2a; 
with eigenvalue huj, where lu = y^ujQ + uif + lu"^, 
|<) = K|0) - a;,|l))/vWTRF and = 

Kc^o*|0)+^aC^r|l)-(koP + kin|a))/(u;Vl^oP + l^iP) 
with eigenvalue 0, and = (-w|e) -I- wJlO) + 

w*|a))/V2i^ with eigenvalue —hw. The WZ phase 
gates are based on ?7°°, in terms of our notation. 
Using the instantaneous eigenstates and eigenval- 
ues, the non-adiabatic correction is obtained as 

En=-i,iJ:^.=o..0 '^AO)Uf{TMm, whose order of 
magnitude is of 1/loT. The two-bit gate proposed there 
is only a Berry phase gate under the Hamiltonian [l^ 

H,, - ^h|f^i|V/iVti + |f^a|V|:at], where a% ^ 
e"^" |e)jj(/z| -I- h.c, (j)i — (j)a = 0/2, using the notations 
therein. The eigenstates are 0^ = (— Ifiipe"*"^!!!) + 
\na\^\aa) + + |r!a|4|ee))/V 2(|»i|4 + 

with eigenvalue ^ y^\ri^]^T\n^ , 

+ \ni\\'^\aa))/./\rkfTW with eigen- 



v alue 0, and = (-jnijV^^l) -I- l^aplaa) - 

V^TI7F^P^|ee>)/V2p7P+FU^ with eigenvalue 

— ^ y^r^ipH-I^XiF- It was proposed to use \(fP) to 
implement the phase gate. The non-adiabatic correction 
is of the order of h/T{rj^ /6)^/\^hJ^T]n^. 

Similar proposals were also made in Josephson junc- 
tion charge qubits 0, For the Hamiltonian 
used in 16], there are an eigenstate with eigenvalue 
\JK^ + I Jip + I J2P, two degenerate eigenstates with 
eigenvalue /i, two degenerate eigenstates with eigenvalue 

— h, which are used to implement the WZ phase gate, and 
one ground state with eigenvalue \Jh? -f | Ji p + | J2 p , 
where h — Ec{l — 2noff)/2. Thus the non-adiabatic 
correction is of the order of h/{^h'^ + | Jip -I- | J2P — 
h)T. In the two-bit implementation, the eigenvalues are 
-V|Jb|2 + (2/i)2, -2h, 0, 2h, V|Jb|2 + (2/i)2. The eigen- 
states with eigenvalue —2h are used as the qubit states. 
The non-adiabatic correction is of the order of h/AT, 
where A is the smaller one of ■\/| Jbp -I- {2h)'^ — 2h and 
2h. Suppose the order of magnitude of Josephson energy 
is J. Then if 1 — 2no// is close to 1, the energy gap is 
of the order of J"^ / Ec in the single-bit gate, and is of the 
order of J^/iE,, in the two-bit gate. Since E^ > J, the 
energy gap is smaller than J. Hence compared with the 
case of [13, the adiabatic condition is harder to meet, 
i.e. the non-adiabatic correction is larger. On the other 
hand, if 1 — 2noff is tuned to be very small, then the 
energy gap for the cases of both and are of the 
order of the Josephson energy. 

In the one-bit gate in the energy eigenvalues 

are 5Ec + [SEc)^ + 2J'^ , 0, which is with twofold de- 
generacy and is used to implement the WZ gate, and 
6Ec — \J {5EcY + J2, where 5Ec is some charging en- 
ergy difference, = \ Jl\^ + \ Jm\^ + \Jr\^ , using the 
notations there. Thus the non-adiabatic correction is 
of the order of h/AT, where A is the smaller one of 
\SEc + VT^^r?TJ2| and \SEc ~ ^/{6EcFTJ^\. Hence 
A is of the order of SEc if SEc > J, and is of the 
order of J if SEc < J- For the two-bit gate, the 

three energy eigenvalues are \J\Jx\^ + |^m^P/2, 0, and 

— \/| JjfP + I j|^''p/2, where the parameters are as de- 
fined there. He nce the non-ad iabatic correction is of the 

order of 2h/^\JxV + VJ^iV^. For both one-bit and 
two-bit gates, the energy gap is at most of the order of 
Josephson energy. Therefore, though WZ gate has more 
advantages over Berry gate, the non-adiabatic error for 
|Ti[l3| is larger than that for [ll| . 

A quantum computing process in a gate array con- 
sists of many gate operations on a large number of 
qubits, hence a complete estimation of error must 
include its scaling with the number of gate opera- 
tions. Suppose from time to T, Mi adiabatic gates, 
denoted as U^{t), (j = 1, - ■ ■ ,Mi), are in parallel 
operation, each on a small number of (say, one or 
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two) qubits. For < t < T, the entire quantum 
computer evolves as = U^{t) ■ ■ ■U^^{t)\'^{0)) = 

En-.M,.. U^mlM ■ ■ ■ U''Hmfi\ (0))|C>, where 
{(pl.) is a basis state of the qubits acted by the j-th 
gate, r denotes the rest qubits, which are not oper- 
ated by any gate during this period. We know that 
U'{m,m = Um^i^m + ^'it), where U^{t) repre- 

sents the adiabatic hmit of j-th gate, while e^ (t) < 1 is 
its non-adiabatic correction. The state of the quantum 
computer at T is then {-^{T)) = |*o(r)) + <t{0,T), 
where 1*0^) = (T) ■ ■ ■ {T)\^ (0)) , 

a{0,T) « is the first-order error of the 

entire quantum computer accumulated from time 
to T. Afterwards, during T < t < 2T, the quantum 
computer is operated in parallel by M2 gates, labelled as 
Ml -|- 1, • • • , M1 + M2, to which the qubits are allocated in 
a way usually different from the period < i < T. Using 
a derivation similar to the above, it can be obtained that 
|*(2T)) = C/o^i+i(T)---[/o^^i+*^^(T)|*(r)) +c7(r,2T), 
where a{T,2T) « EjIm^+'i e^' (r) is the er- 
ror of the entire quantum computer accumu- 
lated from T to 2T. Therefore |*(2r)) = 
jjm+M, . . . u^i'+\T) ■ ■ ■ U^'^' (T) • • • C/oi|*(0)) + 
cr(0,2T), where cr(0,2T) w cr(0,r) -I- a{T,2T) « 
Y^m+M2^jfj.^ is the total error at 2T. Therefore, 
for a quantum computing process consisting of many 
gate operations, no matter how they are arranged in 
space and time, the total non-adiabatic error is, to the 
first-order, just the sum of the errors of all these gates. 
Suppose for each adiabatic gate, the time duration 

< T, and the minimum energy gap with other eigenstates 

< A. Thus the lower bound of the non-adiabatic error 
for each gate is ||e(r)|| ^ %/ I^T. Hence the lower bound 
of the total error is cr ~ Me{T), where M is is total 
number of gate operations. For a quantum computing 
process to make sense, it is constrained that \\a\\ < 1. 
Therefore, M < l/||e(r)|| = AT/h. In Shor's algorithm, 
to factor a number N, M ~ 300(logio Af)^ There- 
fore, adiabatic quantum computing can at most factor 

TV « 10(l/300ll.(T)||)i/3_ of ^jjg o^dg^ 

N « 10. 

Let us switch to the adiabatic quantum algorithm 'sj, 
which is based on adiabatically varying the Hamiltonian 
from a beginning Hamiltonian Hh at t = to a final 
one Hp at t — T. Under the adiabatic limit, if the sys- 
tem starts with the ground state of Hf,, it ends up as 
the ground state of Hp, which gives the solution to an 
optimization problem. For a finite varying rate of the 
Hamiltonian, to the first order, the more accurate state 
is as given in Eq. Q. Hence to the first order approxi- 
mation, the non-adiabatic correction at time T is of the 
order of 1/ApT, where Ap is the energy gap of Hp, which 
is independent of the specific path in which Hh is evolved 
to Hp. 

According to (PJ), as far the perturbative approach is 



valid, i.e. \{(l)oi.t)\dt(f>n{t))/{En - Eo)\ < 1, the mea- 
surement shows that one of the eigenstates appears with 
probability clearly the largest. Then one can know that 
this state corresponds to the ground state and thus the 
solution to the problem. This is consistent with [^. 

To summarize, we developed a general, perturba- 
tive, formulation of the adiabatic quantum computing 
schemes, which perturbatively describe the accurate evo- 
lution of the state. It leads to a deeper understanding of 
related issues. The formalism is applied to analyze both 
the the adiabatic geometric quantum computation and 
the adiabatic quantum algorithm. The order of magni- 
tude of the first-order non-adiabatic error is the inverse 
of the executing time times the minimum gap with other 
eigenstates. Several proposed physical implementations 
of the former are considered from this point of view. 
Different proposals based on charging Josephson junc- 
tions are compared. We also consider an entire quantum 
computing process consisting of many adiabatic gates, 
obtaining the lower bound of the non-adiabatic error, 
as an interesting constraint on the power of the quan- 
tum computation based on adiabatic geometric gates. 
One needs to enlarge the energy gap in order to reduce 
the non-adiabatic error and thus improve the computa- 
tional power. For the adiabatic quantum algorithm, it 
is noted that it can be realized as far as the perturba- 
tive approach, rather than the rigorous adiabatic limit, 
is valid, hence the computational time may be appropri- 
ately shortened. 

Y.S. was supported by a Humboldt Fellowship when 
visiting Germany. 
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